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ABSTRACT 

The paucity of reliable achromatic breaks in Gamma-Ray Burst afterglow light curves 
motivates independent measurements of the jet aperture. Serendipitous searches of af- 
terglows, especially at radio wavelengths, have long been the classic alternative. These 
survey data have been interpreted assuming a uniformly emitting jet with sharp edges 
("top- hat" jet), in which case the ratio of weakly relativistically beamed afterglows to 
GRBs scales with the jet solid angle. In this paper, we consider, instead, a very wide 
outflow with a luminosity that decreases across the emitting surface. In particular, we 
adopt the universal structured jet (USJ) model, that is an alternative to the top-hat 
model for the structure of the jet. However, the interpretation of the survey data is 
very different: in the USJ model we only observe the emission within the jet aperture 
and the observed ratio of prompt emission rate to afterglow rate should solely depend 
on selection effects. We compute the number and rate of afterglows expected in all-sky 
snapshot observations as a function of the survey sensitivity. We find that the current 
(negative) results for OA searches are in agreement with our expectations. In radio 
and X-ray bands this was mainly due to the low sensitivity of the surveys, while in the 
optical band the sky-coverage was not sufficient. In general we find that X-ray surveys 
are poor tools for OA searches, if the jet is structured. On the other hand, the FIRST 
radio survey and future instruments like the Allen Telescope Array (in the radio band) 
and especially GAIA, Pan-Starrs and LSST (in the optical band) will have chances to 
detect afterglows. 



1 INTRODUCTION 



Surveys for transient sources may detect Gamma- Ray Burst 
(GRB) afterglows. In this paper, we call an "orphan" af- 
terglow any afterglow associated with such serendipitous 
searches, as opposed to "triggered" GRB afterglows, local- 
ized through the preceding prompt 7-ray emission. Rhodes 
(1997) suggested that these surveys could be used to put 
constraints on the geometrical beaming angle of the GRB 
jets in the "top-hat" (thereafter TH) model. In this model, 
GRBs are assumed to be uniformly emitting within a cone 
of angle Sjct with sharp edges, where the luminosity drops 
suddenly to an undetectable level. 

The suggestion by Rhodes (1997) is based on the fact 
that the prompt 7-ray emission is relativistically beamed 
within an angle 9jct + 1/r, where l/F << 9jct- Thus, if the 
line of sight lies outside this angle, the GRB is unlikely to be 
detected. However, as the outflow slows down in the after- 
glow phase, the visible region increases to eventually encom- 
pass the observer's line of sight. At late times, F ~ 1 and 
the emission is roughly isotropic. This behaviour would sug- 



gest that more long-wavelength transients than 7-ray ones 
may be expected. Detections of transients at different wave- 
lengths may thus be used to constrain the beaming factor 
b oc . A measure of this quantity would be of great im- 
portance as it would allow one to calibrate Sjet and then 
estimate the true GRB rates (oc 6) and energetics (oc b~^). 

Searches for OAs have been performed at various wave- 
lengths, but none of them has yielded a firm detection. Sev- 
eral authors have used observations in the radio band to con- 
strain the GRB rate and energetics (e.g. Perna & Loeb 1998; 
Woods & Loeb 1999; Paczyski 2001; Levinson 2002; Gal- 
Yam et al. 2006), since late time (i.e. nearly isotropic) af- 
terglow emission peaks in this band. The simplest and most 
common assumption is that the predicted number of OAs in 
a snapshot observation is proportional to the beaming factor. 
Perna & Loeb (1998) used the lack of detections to set an 
upper limit of 6 < 10"^. More recently, Levinson et al. (2002) 
and Gal- Yam et al. (2006) compared their data with a more 
detailed model for radio afterglows. They showed that the 
number of expected OAs in a flux limited survey is inversely 
proportional to b and they place a lower limit of 6 ^ 60. 
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X-ray survey data have been used to search for OAs 
by Grindlay (1999) and by Greiner et al. (2000), while op- 
tical searches have been more numerous (Schaefer et al. 
2002; Vanden Berk et al. 2002; Becker et al. 2004; Rykoff 
et al. 2005; Rau et al. 2006, Malacrino et al. 2007). Using 
shorter wavelengths than radio to constrain the beaming 
factor necessarily requires a careful comparison with theo- 
retical predictions (e.g. Totani & Panaitescu 2002; Nakar, 
Piran & Granot 2002), since the emission is likely to be 
still relativistically beamed when OAs are detected in those 
bands. Malacrino et al. (2007, 2007b) put the tightest op- 
tical constraints so far (see their fig. 3). Their non-detection 
resulted in an upper limit for the number of OAs on the sky 
that is marginally consistent with the predictions of Totani 
& Panaitescu (2002, thereafter TP02) and consistent with 
Nakar et al. (2002, thereafter N02) and Zou et al. (2007, 
thereafter Z07). 

The purpose of this paper is to predict results for 
searches of afterglows in surveys, for a different GRB jet 
structure. Currently, OA surveys generally aim at constrain- 
ing the jet angle. However, if GRB outflows are not geomet- 
rically beamed but they rather have an anisotropic luminos- 
ity distribution, the interpretation of data within the TH 
scenario would be misleading. 

Numerical simulations of collapsing massive stars (e.g. 
MacFadyen & Woosley 1999) show that the jet emerges from 
the star with an energy distribution E{0) and Lorentz factor 
r{6) that vary as a function of the angle 6 from the jet axis. 
It has been shown (Rossi, Lazzati & Rees 2002; Zhang & 
Meszaros 2002) that, if E(e) oc the diversity of after- 
glow light curves can be ascribed to different viewing angles 
within the context of an universal structured jet (US J). In 
the USJ model, the outflow is geometrically wide. In this 
paper, we will postulate that for each GRB two simultane- 
ous and opposite jets with 9jct ~ 90° are produced. In this 
model & = 1. The feature of having emission into Att solid 
angle is attractive since it can explain the lack of OA detec- 
tions, while the non-uniform energy distribution allows one 
to avoid the huge energy requirement, demanded by GRBs 
with an isotropic equivalent energy ^ 10^* ergs (e.g. GRB 
990123). An important consequence of the energy distribu- 
tion law E(d) oc is that it establishes a unique relation 
between the viewing angle and the observed luminosity, once 
a radiation efficiency law with angle is assumed. Thus, un- 
like in the TH model, the observed luminosity function is 
not a free parameter. Consequently, the uncertainties in the 
predicted OA rates in the USJ scenario are smaller than in 
the TH model (see §15). 

In this paper, we use the USJ framework to compute 
the expected number of transients in an all-sky snapshot 
and their rate as a function of the survey sensitivity. The 
aim is pursued by means of Monte Carlo simulations of the 
afterglow properties, as observed by X-ray, optical and ra- 
dio surveys. The procedure is described in § [21 We show 
prospects for OA detections with current and future surveys 
in § [3l This allows us to identify the best survey character- 
istics to increase the chance of detection and single out the 
most promising future missions. A comparison of our results 
with the top-hat predictions is performed in §[4] Finally, we 
discuss and conclude our work in § [S] 

Throughout the paper, we assume a flat Universe with 
Ho = 73 km s"^Mpc"\ D,m = 0.3 and Ha = 0.7. 



2 SIMULATING THE POPULATION OF 
ORPHAN AFTERGLOWS 

We use Monte Carlo methods to simulate the population of 
GRB orphan afterglows in the USJ. GRBs are randomly gen- 
erated on the sky with a probability distribution in redshift 
that traces the star formation rate (SFR) (§[2111 • We use the 
external shock model to compute the afterglow luminosity 
curve (§[2]2]). The probability function for the viewing angle 
6 is given by the fraction of the solid angle associated with 
that angle, P{0) oc sin(S). Our simulation yields for radio, 
optical and X-ray bands the distribution of afterglow fluxes 
and the total number of afterglows on the sky for a snapshot 
observation, together with the average time Tth that an af- 
terglow remains detectable in the sky, as a function of the 
detection threshold. Finally, we compute the OA detection 
rate for any flux limited survey where the observation time 
is much greater than Tih- 

2.1 Formation rate and 7 ray luminosity function 

We assume that the GRB population in the universe is de- 
scribed by a redshift-independent luminosity function, a red- 
shift distribution and a distribution of the spectral parame- 
ters. 

In the USJ, the isotropic equivalent kinetic energy in the 
afterglow phase has the angular dependence 
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where the bright central spine with angular size 0c has a 
maximum kinetic energy Ec — E{0). The expected lumi- 
nosity function follows from Pgi<.b{L) — P{d)^ (Rossi et 
al. 2002), 
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where L [erg s~^] is the isotropic equivalent bolometric peak 
luminosity and Lc is proportional to Ec (see eq. [4]). In equa- 
tion[2l we assume that the 7— ray emission efficiency and the 
ratio of mean luminosity to peak luminosity is independent 
of the angle. 

The GRB comoving rate -Rgrb(2) [yr~^ Mpc"''] is assumed 
to follow the comoving rate Rsn{z) [yT~^ Mpc"''] of Type II 
supernovae -Rgrb(z) — k x Rs-n{z), where k = Rgrb/Rsn 
is a free parameter of the model. We assume that these arise 
from stars with masses above 8 Mq and that the initial mass 
function has a Salpeter form (e.g. Porciani & Madau 2001). 
Daigne, Rossi & Mochkovitch (2006, thereafter DRM06) 
found that the above prescription is dubious for redshift 
greater than 2. However, most of the OAs that are detected 
in a survey are located at lower redshifts, where this assump- 
tion appears to hold. In our analysis, we find, in fact, that 
the mean OA redshift for any reasonable ffux threshold is 
never greater than z = 2. The star formation rate we adopt 
to derive _Rsn (dubbed SFR2) and its comparison with data 
(Hopkins 2004) are shown in fig. 1 of DRM06. This SFR 
saturates beyond z ~ 2 at a level of 0.2 Mq yr~^ Mpc"''. 
Even though this behavior is consistent with data at the 
one sigma level, the flat extrapolation of our SFR after the 
peak remains questionable, since high-z data are plagued by 
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uncertainty on the amount of dust extinction. However, as 
discussed above, we do not expect that uncertainties in the 
high-z behaviour of Rgkb{z) appreciably affect our results. 
The GRB redshift probability function is thus given by 



dV Rgkb{z) 
dz 1 + z 



(3) 



where the comoving volume equals 



dV _ c 4ttDI{z){1 + z)-'^ 
~ Ho (1 + + 

and Di^{z) is the standard luminosity distance. 
The spectral properties of GRBs are described by the distri- 
bution of their peak energy Ep and their low and high-energy 
slopes a and (3. For the peak energy, we assume a log-normal 
distribution with a mean value £p,o and a dispersion 0.3 dex. 
For the slopes, we adopt the observed distribution by Preece 
et al. (2000). 

With these assumptions, the GRB population is en- 
tirely described by four free parameters : Lc and 9c for the 
luminosity function, k for the comoving rate and Ep^ for 
the spectral properties. Following the method described in 
DRM06, we constrained them by fitting simultaneously: i) 
the logN — logP distribution of GRBs (where P [ph cm~^ 
s~^] is the peak flux) detected by the Burst and Transient 
Source experiment (BATSE; Kommers et al. 2000; Stern et 
al. 2000, 2002); ii) the peak energy distribution of bright 
BATSE bursts (Preece et al. 2000); and ui) the HETE2 
fraction of X-ray rich GRBs and X-ray flashes (Sakamoto 
et al. 2005). For a given set of parameters {Lc,9c,k, Epfi), 
a population of ^ 10^ GRBs is randomly generated, with a 
redshift z, a luminosity L and a spectrum characterized by a 
peak energy Ep and a low and a high-energy slope a and (3. 
The four free parameters are then adjusted to minimize the 
obtained when comparing the simulated data with the 
three observations listed above. More details about the pro- 
cedure can be found in DRM06. The results for the best fit 
are log(Lc[erg s'^]) = 53.7 ±0.6 and Oc = 9.2 ± 5.2° for the 
luminosity function, log(k) = —5.99 ± 0.06 for the comov- 
ing rate and log (i5p,o[keV]) — 2.8 ± 0.1 for the spectrum. 
The reduced = 1-53 for 37 degrees of freedom. We note 
here that DRM06, assuming a power-law luminosity func- 
tion, found ^ —1.6 as the best fit value for the slope when 
this is free to vary, as opposed to the slope ~ —2 predicted 
by the US J model (eq. [2]). However, an acceptable is also 
obtained in our case. 



2.2 Physical description of afterglow lightcurves 

The afterglow emission is modeled as synchrotron radiation 
from a relativistic blast-wave propagating in a constant den- 
sity external medium (e.g. Meszaros & Rees 1997). We ig- 
nore the contribution from the reverse shock, although it 
dominates the afterglow emission in the first few tens of 
seconds after the GRB (Sari & Piran 1999) , since OA obser- 
vations occur much later. For the same reason, our results 
are largely independent of the choice of the initial Lorentz 
factor (which we fix at Fq = 300), since the deceleration 
time is unlikely to exceed 10"^ s (e.g. Panaitescu & Kumar 
2000). Similar considerations of timing allow us to neglect 
the modelling of the early afterglow features (flares, plateau, 
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Figure 1. Flux distribution in the X-ray band at 1 keV energy; 
these distributions are intrinsic, i.e. no selection criteria have been 
applied. We choose the minimum flux, so that the arbitrary cut- 
off at 10 yrs for the age of an afterglow in our simulations does 
not affect the shown distributions. Left panel: the distribution 
of fluxes at 12 hours after the trigger. This may be compared 
with fig. 5 of Berger et al. (2005), taking into account that their 
observed distribution suffers from selection eSects. Right panel: 
the flux distribution as it appears in a snapshot observation of 
the sky. 
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Figure 2. As Fig. [T] but for magnitude distributions in the R 
band. Our left panel may be compared with fig. 4 of Berger et al. 
(2005) 



etc., e.g. Nousek et al. 2006), that are not reproduced by 
the standard external shock model. 

The code we use for the calculation is described in Rossi 
et al. (2004). For each lightcurve. La [erg s~'^Hz~^], the in- 
put parameters are: the angular distribution of kinetic en- 
ergy, the viewing angle, the shock parameters, the external 
density and the rest-frame frequency. 

The luminosity function parameters, Lc and 9c, con- 
strained by the prompt emission data (§[2TT} allow us to de- 
termine the kinetic energy distribution within the jet (eq.[l|, 
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Figure 3. As Fig. [T] but for flux distributions in the radio band 
at 5 GHz frequency. Our left panel may be compared with fig. 6 
of Berger et al. (2005) 
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Inspection of current data suggests that the proportional- 
ity constant is of the order of 1 second. The 7— ray effi- 
ciency rj.y that is inferred from modehng data in the stan- 
dard external-internal shock scenario is rather high, be- 
tween ~ 50% and ~90% (Panaitescu & Kumar 2002). This 
uncertainty together with the one a range of imply 
lO'^^ < erg < 2 X 10^^. We adopt = 1.3 x 10'^=' erg. 
Since the core angle is not strongly constrained by prompt 
emission data (see 12. ip . we adopt a value at the low end 
of the one a range, 9c =4°, to account for the observed 
breaks in the lightcurves on day timescales. The external 
shock and density parameters are chosen from the ranges of 
values inferred from afterglow modeling (e.g. Panaitescu & 
Kumar 2001; Panaitescu 20050 . The fraction of energy at 
the shock that goes into accelerated electrons and magnetic 
field is Ee = 0.05 and en = 0.005 respectively; the electrons 
are accelerated into a power-law with exponent p = 2.2. 
The external number density is taken to be n = 1 cm~'^. 
The whole set of parameters {Ec, 9c, te, es, p and n) yield 
results consistent with the observed afterglow flux distribu- 
tions (Berger et al. 2005). Our flux distributions are shown 
in Fig. [1] to Fig. [3] In the right panels, we plot the his- 
togram of fluxes in a given band for all afterglows detected 
in a snapshot observation of the sky. In the left panels, we 
show the flux histogram of the same afterglows evaluated 
at a common observed time. These latter distributions com- 
pare favorably with fig. 4, fig. 5 and fig. 6 of Berger at al. 
(20050. 



^ We are aware that shock and density parameters inferred from 
observations are not universal. They rather vary from burst to 
burst within some ranges. To model the emission, however, what 
is important is that the combination of these parameters, includ- 
ing the fireball kinetic energy, gives fiuxes comparable to observa- 
tions. As mentioned afterwards, we perform such a comparison. 
It indicates that we have chosen a reasonable combination of pa- 
rameters. 

^ We do not attempt a formal quantitative comparison with data, 



We compute afterglow lightcurves I/a(i' x (1 + z),9,t') 
for three observed frequencies, u = 2.42 x 10^^ Hz (1 keV), 
1/ = 4 X 10" Hz (R band) and = 5 x 10^ Hz; 6 redshifts, 
^z^ 20; 10 viewing angles, 0° ^ 9 ^ 90° and 57 comoving 
times t' , spanning 10 years. We arrange those data in the 
form of a matrix that can be easily interpolated in order 
to assign a luminosity La to each simulated burst. The flux 
observed on Earth is calculated from La with F — La(l -I- 
2:)/(47rD£). Examples of lightcurves are shown in Fig. |3] in 
this example the GRB is situated at z = 1 and viewed under 
different angles, in our three bands of observation. 



2.3 Monte Carlo code for "Orphan" Afterglows 

OAs are generated with a flat probability distribution in 
age, ta, (i.e. the time lag as observed on earth between the 
GRB explosion and the snapshot observation epoch), in an 
interval of 10 years, at a rate of Rohe — 2195 yr~^ 0- For 
the sensitivities of interest in this paper, afterglows of ten 
years or older are undetectable. Their mean duration above 
threshold is indeed shorter than our total simulated time 
(see Figs. [5] [6] and [7)1. Therefore, the arbitrary cut off at ten 
years does not influence our results. 

We generate each GRB redshift, z, according to the 
probability distribution discussed in § 12.11 and its viewing 
angle 9 according to 



P(< 6*) = (1 - cose). 



(5) 



with 0° ^ 9 ^ 90°. From the matrix of the afterglow lumi- 
nosities interpolated at ta/ {1 + z), 9 and 1/ {1 + z),we calcu- 
late the OA flux and compare it with a given flux threshold, 
Fth- We can, thus, compute the expected number of after- 
glows above Fth in a given band ly: Nsnap{> Fth, i')- This is 
the number that would be seen in a snapshot observation of 
the entire sky. 

Model predictions for a flux limited survey require the 
calculations of two other quantities. First, the mean of the 
time interval tth spent by an afterglow above the flux limit. 



{logio ith) = = logio Lth, 



(6) 



where we actually compute the geometric mean of tth to 
avoid being biased by extreme values. Second, the mean rate 
at which afterglows appear in the sky over the survey flux 
threshold. 



Roa. — Nsi 



'■th 



iVs, 



(7) 



We performed 160 simulations of flux limited all sky 
snapshots in radio and optical bands and 400 in X-rays, 
where we needed more statistic. We then computed the av- 
erage A^snap, Irate and Tth . The results are shown in Figs O [5] 
and [7] (thick lines). They can be used to estimate the num- 
ber of OAs expected in a given survey. If the observing time 
of a survey Tobs is shorter than the time Tth for which an 



since this would require us to take into account selection effects 
in diff'erent bands that are difficult to quantify. 

This expected rate of GRBs observed from the Earth is given 
by integrating -Rgrb/(1 + z) = feiJsN/(l + ^) over the whole 
volume of the universe. 
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afterglow is detectable, then we can consider it a snapshot 
observation and the total expected number of detected OAs 
is 

iVoa(> F,l,, V) = iV.nap(> Fth , v) , (8) 

47r 

where f2oba is the solid-angle of the sky covered by the snap- 
shot. Vice versa, the total expected number of OAs is com- 
puted as 

Afoa(> Fth, V) = 7?oa(Fth, V) Tobs ■ (9) 



3 RESULTS 

Our results show that, as expected, the probability of detec- 
tion and the mean duration above threshold increases with 
the survey sensitivity (Figs. [5l [6] and [7)1 . A chance of detec- 
tion in an all sky snapshot (A'snap ^ 10) requires: a flux limit 
of vFv 10"^'* erg cm~^ s~^ at 1 keV ; a limit magnitude 
of 7? ^ 23 in R band and a flux density threshold of ^ 1 
mjy at 5 GHz. 

3.1 Specific surveys predictions 

In this section, we provide illustrative predictions for var- 
ious current and planned surveys. We use here the results 
shown in Figs [5] |6] and [7] and eq.[8]and eq.|9l We note that 
the following numbers should be taken as upper limits when 
compared with observations, since we did not consider de- 
tection limitations (e.g. host galaxy, dust absorption etc..) 
other than instrumental. 

We also note that a detailed comparison between obser- 
vation and theory requires knowledge of the specific survey 
strategies. However, generic conclusions can be drawn from 
the following examples. 

3.2 X-rays 

The two most sensitive X-ray surveys (~ 1 x 10^^^ erg 
s~^ cm^^ in the 0.5-2 keV band) performed by Chandra 
and XMM-Newton are the 2 Msec Chandra Deep Field- 
North and the 0.8 Msec XMM-Newton Lockman Hole field. 
They cover respectively an area of 0.13 and 0.43 deg^. Since 
Tth — 52 days, these surveys are equivalent to two snapshot 
observations. The small coverage of the sky results in a very 
small detection probability of a few 10~*. Larger surveys, 
as the XMM-Newton Bright Serendipitous Source Sample 
(Delia Ceca et al. 2004) or The Chandra Multi-wavelength 
Project (ChaMP; e.g. Kim et al. 2007), are less sensitive 
(I'Fth ~ 10~^* erg s~^ cm~^) and the increase in area is not 
sufficient to yield more than A'^oa ~ 10~^. 

The ROSAT AU Sky Survey (RASS; e.g. Voges et al. 
1999) covers the full sky. The RASS exposure is 76435 
deg^ days, with a sensitivity of 10~^^ erg s~^ cm~^. There- 
fore the survey is equivalent to an all sky observation with 
Jobs = 1.85 days. At this fiux threshold, OAs are fast X- 
ray transients, lasting for Tth — 0.5 days and Roa. — 0.1 
day~^. Therefore, despite the large coverage of the sky, we 
predict that the survey should have found only A^oa — 0.2. 



Greiner et al. (2002) found 23 OA candidates. After spec- 
troscopic follow-up, however, they concluded that most, if 
not all events, are stellar flares. This is in agreement with 
our predictions. 

The prospects for detection are not exciting even for 
the future mission eROSITA (extended ROentgen Survey 
with an Imaging Telescope Array). It will perform the first 
imaging all-sky survey up to 10 keV with a sensitivity of 
5.7 X 10"" ergs"^ cm"^ in the 0.5-2 keV band. We predict 
~ 2.3 OAs. 

These dispiriting results can be understood from Fig. [1] 
and Fig. [S] the flux limit should be ^ 10"^"* erg s~^ cm~^ 
in order to have a good chance of detection in one all sky 
snapshot. To achieve such a sensitivity, a long exposure time 
is needed (e.g. ~ 50 ksec for Chandra). Given the small 
field of view of the current instruments, a full sky scan is 
unfeasible. We conclude that X-ray surveys are a poor tool 
for OA searches if jets are described by the USJ model. 

3.3 Optical 

The most recent optical OA searches are the ones by Ran et 
al. (2006) and by Malacrino et al. (2007). Rau et al. observed 
12 deg^ of sky for 25 nights, separated by one or two nights. 
They used the MPI/ESO Telescope at La Silla, reaching 
R=23. Since at that sensitivity Tth — 22 days and Trate — 2 
days, we compute the expected OAs through an average 
rate of -Roa ~ 6.3 day~^. We find A'oa = 0.02, considering 
an actual observing time of 12.5 days. Malacrino et al. used 
images from the CFHTLS very wide survey . They searched 
490 deg^ down to R=22.5. They observed 25-30 deg^ every 
month over a period of 2-3 nights (Malacrino et al. 2006). 
Since Tth — 17 days, we can consider their survey a snapshot 
observation and we predict A'oa = 0.1. Recently, Malacrino 
et al. (2007b) have rejected the only candidate they had 
identified in their first work. 

The most recent Sloan Digital Sky Survey (SDSS) data 
release (Adelman-McCarthy et al. 2007) includes imaging 
of 9583 deg^, with a R magnitude limit of 22.2. We expect 
Aoa — 1.5. If we take a more conservative magnitude limit 
of 19 to account for the need for spectroscopic identification, 
the detection probability drops to ~ 0.05. 

The previous meagre results are due to the fact that a 
very large detection area is essential for these sensitivities 
(fig. fright panel and fig. O. Only with the fiux limit of 
the Subaru Prime Focus Camera (R~ 26 in 10 minutes of 
exposure), we could restrict ourselves to 5 % of the sky and 
get a snapshot with A^oa — 10. This, however, would require 
a total observing time of Tobs — 57 days. 

Future larger surveys include GAIA (Parryman et al. 
2000) and the Panoramic Survey Telescope & Rapid Re- 
sponse System (Pan-Starrs; e.g. Kaiser et al. 2002). The first 
is an all sky survey with a magnitude limit of i? ~ 20. It 
will observe each part of the sky 60 times separated by one 
month (Lattanzi et al. 2000). At this sensitivity, an OA stays 
in the sky on average for ~ 3.8 days; thus GAIA will perform 
60 independent snapshots of the sky, with a prediction of 
A'oa — 44. Pan-STARRS is expected to scan three-quarters 
of the entire sky in about a week, down to an apparent mag- 
nitude of 24. Since Tth — 43 days, this can be considered 
a snapshot observation and we expect ~ 23 OAs. Finally, 
the Large Synoptic Survey Telescope (LSST) is planned to 
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Figure 4. Afterglow lightcurves for a GRB at z=l in X-rays (upper left panel), optical (upper right panel) and radio (lower left panel) 
bands. For each band, the different curves correspond to different viewing angles. From top to bottom, 9 = 0°,4°,8.7°, 19° , 41.2°, 90° . 



cover 10,000 square degrees every three nights down to a 
depth of R~ 24.5 (Ivezic et al. 2008). This would yield ~ 13 
afterglows every three nights. Both Pan-STARR and LSST 
will be able to repeatedly observe an afterglow source, thus 
monitoring its variability and enhancing the chances of iden- 
tification. 

Thus, future optical surveys could be powerful tools for 
OA searches. 

3.4 Radio 

Levinson et al. (2002) searched for transients by compar- 
ing the FIRST and NVSS (NRAO VLA Sky Survey) radio 
catalogues and found 9 candidates. Gal- Yam et al. (2005) 
rejected all candidates by means of follow-up radio and opti- 
cal observations and placed an upper limit (95% confidence) 
of 65 radio transients for the entire sky above 6 mjy at 

1.5 GHz. This may be translated to a sensitivity thresh- 
old of 3.3 mJy at 5 GHz, using a typical late time spectral 
shape in radio of _F oc u^"'^ (TP02). We predict ~ 1.4 radio 
afterglows. 

Recently, Bower et al.(2007) published an archival sur- 
vey with data from the Very Large Array, spanning 22 years. 



For an effective area of 10 deg^, we get a rate of ~ 4 x 10~^ 
yr"'^ for OAs brighter than 370 /iJy. This rate is too low to 
account for the 10 detected transients. In addition the ob- 
served transient duration of approximately a week suggests 
that those sources are not indeed afterglows, which are ex- 
pected to last above that threshold for approximately half a 
year. 

Fig. [3] and Fig. [7] and the above examples show that the 
OA search would greatly benefit from lowering the survey 
sensitivity below 1 mJy, with an area of ^ 10,000deg^. 

FIRST (Faint Images of the Radio Sky at Twenty-cm; 
Becker, White & Helfand 1995) has covered over 10" deg^ of 
the North Galactic Cap. The survey area has been chosen 
to coincide with that of the SDSS. The sensitivity is Fth 1 
mJy at 1.4 GHz. This may be translated into a flux limit of 
0.5 mJy at 5 GHz. If we consider, as for the SDSS, an area 
of 9583 deg^, we expect ~ 7 OAs. 

The plan is for the Allen Telescope Array (ATAj3 to 
observe 10** deg^ at mJy sensitivity and later to go as deep 



" See e.g. |http://ral. berkeley.edu/ata/science/ 1 
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Figure 5. Expected number of afterglows in an all-sky snapshot observations in tlie X-ray band, as a function of the flux threshold 
(solid thick line). The thin solid lines are the la contours. We also plot Tt^ (thick dashed line) and Trate (thick dot-dashed line) in days. 



as ~ 0.1 mjy at 5GHz. The expected number of OAs would 
then rise from a few to 50. 



4 COMPARISON WITH THE "TOP-HAT" 
MODEL 

In fig. IHl we compare our results for A'^snap (solid line) with 
published predictions for the TH model in R-band (dashed 
lines). From bottom to top, we plot the "standard" model 
by Z07; the "preferred" and the "optimistic" models by N02 
and the model by TP02. The curves from these works have 
similar slopes but they vary in normalisation by several or- 
ders of magnitude. The main difference is the assumed af- 
terglow luminosity function. 

This is the result of different choices for the opening 
angle and total jet energy distributions and the afterglow 
radiation efficiency. Z07 assume constant total (beamed cor- 
rected) energy Etot = 10^^ erg and a power-law distribu- 
tion of opening angles P(6'jct) oc 6^1- N02 assume constant 
peak flux at 6ohs = ^jot and a fixed average opening an- 
gle 6 jet =0.1 rad ("canonical" model) and 9jat ~ 0.05 rad 



("optimistic" model). Finally, TP02 assume that the whole 
GRB population is represented by 10 well-studied events. 
We also note that within each model there is an order of 
magnitude uncertainty. For example, Zou et al. results dif- 
fer by an order of magnitude between Etot = 5 x 10^" erg 
and Etot = 5 x 10^^ erg (their fig. 4). 

These uncertainties in the TH model are due to the 
fact that neither the opening angle distribution nor the to- 
tal jet energy distribution is uniquely defined by the model 
or unbiasely determined by observations (e.g. through the 
observed luminosity function). The US J model, instead, has 
more predictive power, since it has less degrees of freedom. 

From the comparison in fig. |8l one concludes that our 
model predictions fall roughly in the middle of the range of 
published results for the TH model. The differences in the 
predictions between our model and any individual TH model 
are large enough to be significant, but there is a sufficiently 
wide class of plausible TH models to encompass almost any 
luminosity function of OA that might be observed in the fu- 
ture. In principle, therefore, OA observations could rule out 
the USJ model (noting of course that there is a similar level 
of uncertainty to our predictions as individual TH models) , 
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Figure 6. The same as Fig. [S] but for the R band. 



but the same cannot be said for the TH model. The TH 
model could be tested by a combination of orphan afterglow 
observations and independent constraints on the parameters 
of the model. 



ous bands depends on the survey strategy, the ratio should 
tend to unity (6 — 1) if events are detectable at arbitrarily 
low fluxes. 



5 DISCUSSION AND CONCLUSIONS 

The realization that GRBs may be jetted has triggered stud- 
ies of orphan afterglows. Early studies assumed that the 
prompt GRB emission comes from a sharp jet, in which case 
the ratio of afterglows to that of GRBs yields a constraint 
on the GRB beaming fraction (or equivalently the opening 
angle of the GRB jet). This parameter is of importance for 
a proper assessment of GRB rates and energetics. 

The structured jet model has offered an equivalent ex- 
planation for afterglow phenomenology. However, its inter- 
pretation of the lightcurve breaks is different: they would 
arise from viewing angle effects and not from geometri- 
cal coUimation. In fact, if GRB jets are indeed structured, 
most, if not all, afterglows should be generally preceded by a 
prompt emission pointing towards the observer. Therefore, 
even if in practice the relative number of detections in vari- 



In this paper, we have investigated the detection 
prospects of afterglows for flux limited surveys, in the USJ 
framework. We conclude that large sky coverage is essential 
in all bands. In addition. X-ray and radio instruments should 
push their flux limit below 10"^"* erg s~^ cm~^ (at 1 keV) 
and ~ 1 mjy (at 5 GHz) respectively. Current and planned 
X-ray surveys are thus not suited for OAs searches, if the 
jet is structured. The FIRST and the future ATA projects 
could be successful in detecting radio OAs. The potential is 
even better for future optical all-sky surveys, such as GAIA 
and Pan-Starrs. We also note that it would be worthwhile to 
exploit the great sensitivity of the Suprime-Cam, for which 
5-10% of the whole sky is sufficient for positive detections. 
Certainly, a combination of X-ray, optical and radio observa- 
tions would yield the most of information and will help un- 
derstanding whether the USJ describes correctly the struc- 
ture of the jet. 
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